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We report a Dalitz-plot analysis of the charmless hadronic decays of charged B mesons to the final 
state i^ ± 7r T 7r ± . Using a sample of (383.2 ± 4.2) x 10 6 BB pairs collected by the BABAR detector, we 
measure CP-averaged branching fractions and direct CP asymmetries for intermediate resonant and 
nonresonant contributions. We find evidence for direct CP violation in the decay B ± — > p°(770)K ± , 
with a CP- violation parameter A C p = (+44 ± 10 ± 4 t? 3 )% . 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 
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I. INTRODUCTION 

The Kobayashi-Maskawa mechanism of CP nonconser- 
vation [l[ has, in recent years, been confirmed through 
observations of direct CP violation in the kaon sys- 
tem [E Q and of both mixing-induced 0, Q and di- 
rect [6J CP violation in the B meson system. How- 
ever, it is striking that despite an enormous experimen- 
tal effort [7J, l8(, no observation of CP violation in the 
decay of any charged particle has yet been made. In 
some cases, for example, in B ± -> 7 r if ± [EEll, the ab- 
sence of an asymmetry is difficult to understand theoret- 
ically 0, [IE EE HEED!. The search for such effects 
is therefore a priority for studies of the weak interaction. 

Several theoretical investigations have suggested 
that large CP asymmetries may be observed in 
B ± -» p°(770)K ± decays [H EE EE M HE HI HE 
[241 . |25j |. Because of the large width of the p° meson, 
this channel must be studied using a Dalitz-plot analy- 
sis of the K^n^n^ final state. Previous measurements 
found A CP (B ± -> p°(770)K ± ) = (+32±13±6±®)% [H 
and A CP {B ± -> p°{770)K ± ) = (+30 ± 11 ± 2 ±f)% H3 
(where the errors are statistical, systematic, and model- 
dependent, respectively), indicating evidence for CP vi- 
olation in this decay, and leaving a strong need for more 
precise and conclusive studies. 

Several other considerations motivate a precise analy- 
sis of _B ± — » if ± 7r T 7r ± decays. The CP asymmetries in 
the decays B+ -> K*°(892)tt+, B+ -> if *°(1430)7r+, 
and B + — > i\"2°(1430)7r + are predicted to be negligi- 
ble [IE H3| compared to the current precision, since 
these are mediated by b — > s loop (penguin) transitions 
only, with no b — > u tree component. A significant CP- 
violation effect in any of these channels would there- 
fore provide a signature of new physics. On the other 
hand, decays such as B + — * p°K + may have compara- 
ble contributions from both tree and penguin amplitudes, 
and their interference is sensitive to their relative weak 
(CP violating) phase difference 7. Results from ampli- 
tude analyses of B^ — > if ± 7r =F 7r ± and other B — ► Kitit 
decays [IE HE H(| can be combined to constrain the 
hadronic parameters and allow a relatively clean deter- 
mination of 7 [33, [33, [33|, 113, [3f| |36| . 

Understanding the dynamics of the B^ — > if ± 7r T 7r ± 
Dalitz plot is also a priority. Previous studies have left 
unresolved the nature of a structure peaking at invari- 
ant mass around 1300 MeV/c 2 in the spectrum 
(denoted /x(1300) in this work). Similar enhancements 
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have also been seen at low K + K invariant masses in 
B -> KKK HI HE HI and B ± -> if ± J ftT = F 7r ± go| de- 
cays. The origin of these structures has aroused consid- 
erable interest among theorists HE HE EE EE El , as it 
is of great importance in the understanding of low energy 
spectroscopy [46j |. 

In this paper we present results from an amplitude 
analysis of B ± — > if ± 7r =F 7r ± decays based on a 347.5 fb _1 
data sample containing (383.2 ± 4.2) x 10 6 BB pairs 
(N B g). Compared to our previous publication (26[, we 
have increased the data sample by 70%, included several 
improvements in reconstruction algorithms that enhance 
the signal efficiency, made several modifications to the 
analysis to increase the sensitivity to direct CP- violating 
effects (for example, by including more discriminating 
variables in the maximum likelihood fit), and improved 
our model of the Dalitz-plot structure. Moreover, we 
have developed a novel parametrization of the coefficients 
used in the fit that ensures good statistical behavior of 
the fitted parameters. 

The data were collected with the BaBar detector [47j 
at the SLAC PEP-II asymmetric-energy e + e~ storage 
rings [48| operating at the T(AS) resonance with center- 
of-mass (CM) energy of y/s = 10.58 GeV. An additional 
total integrated luminosity of 36.6 fb _1 was recorded 
40MeV below the T(4S) resonance ("off-peak" data) and 
was used to study backgrounds from continuum produc- 
tion. 



II. AMPLITUDE ANALYSIS FORMALISM 

A number of intermediate states contribute to the de- 
cay B^ — ► K^n^n^ . Their individual contributions are 
obtained from a maximum likelihood fit of the distribu- 
tion of events in the Dalitz plot formed from the two 
variables m 2 Kv = m 2 (±nT and m 2 ^ = m^ ±jT . The total 
signal amplitudes for B + and B~ decays are given, in 
the isobar formalism (see for example [IE HEHU), by: 

A = A{m 2 K ^ml n ) = ^^(m^.mjj (1) 

i 

A = A(m 2 Kn ,ml^) = ^ CjFj(m 2 K7T , rojj, . (2) 
3 

The complex coefficient for a given decay mode j is Cj 
and is measured relative to one of the contributing chan- 
nels (if* (892) in this analysis). The Cj contain all the 
weak phase dependence and so Fj = Fj. The distribu- 
tions Fj describe the dynamics of the decay amplitudes 
and are the product of an invariant mass term (Rj ) , two 
Blatt-Weisskopf barrier form factors (Xj(z)) [52j, and an 
angular function (Tj): 

Fj = Rj x Xj(p*) x X. ; (q) x Tj , (3) 

where J is the spin of the resonance, q is the momentum 
of either daughter in the rest frame of the resonance, and 



5 



p* is the momentum of the bachelor particle in the rest 
frame of the B. The Fj are normalized over the entire 
Dalitz plot: 



\Fj(rn? Kw ,ml w )\ dm 2 K7r dm 



2 

7T7T 



1 . 



(4) 



The Blatt-Weisskopf barrier form factors are given by: 

X.j=v(z) = 1, (5) 
X J= i{z) = \A/(1 + (^bw) 2 ), (6) 
X J=2 {z) = Vl/((^^Bw) 4 + 3(zrBw) 2 +9), (7) 

where rew, the meson radius parameter, is taken to be 
(4.0 ± 1.0) (GeV/c)" 1 [E3. 

For most resonances in this analysis the Rj are taken 
to be relativistic Breit-Wigner line shapes: 



Rj(m) 



(tjIq — to 2 ) — imoTim) 



(8) 



where mo is the nominal mass of the resonance and T(to) 
is the mass-dependent width. In the general case of a 
spin- J resonance, the latter can be expressed as 



r(m) 



r„ ^ 



, 2J+1 v2 / \ 

>1 \ nj}0\ Xj(q) 
q J \ m J Xj(q )' 



(9) 



The symbol Tq denotes the nominal width of the reso- 
nance. The values of too and To are obtained from stan- 
dard tables Q. The symbol qo denotes the value of q 
when to = too- 

For the /o(980) line shape the Flatte form [54| is used. 
In this case the mass-dependent width is given by the 
sum of the widths in the 7T7t and KK systems: 



where 



r(m) = r 7rir (m) + T KK (m), 



r 7rTr (m) = gA-^Jl- 4ml /m 2 + 



(10) 



(11) 



9K 2 A/1 



1 — 4m 2 + lm 2 



4m 2 K±l 



/m 2 - 
i y/l - Am 2 K0 /m 2 



(12) 



The fractional coefficients arise from isospin conservation 
and g n and gK are coupling constants for which we take 
the values: 

. 9rr = (0.165 ± 0.010 ± 0.015) GeV/c 2 , (13) 
g K = (4.21 ± 0.25 ± 0.21) x 5jr , 

from results obtained by the BES experiment [55j . 



The + component of the Ktt spectrum is not well un- 
derstood [H, |57|; we dub this component (i^7r)Q° and 
use the LASS parametrization (5(| which consists of the 
i^o°(1430) resonance together with an effective range 
nonresonant component: 



Rj(m) 



q cot 8 b — iq 



(14) 



TOqTi 



mo 



(ml - mj <1T ) - iTO F c 



mK-K qo 



where cot 8b — + \rq. We have used the following 
values for the scattering len gth and effective range pa- 
rameters of this distribution [261 ]: 



a = (2.07 ± 0.10) (GeV/cT 1 
r = (3.32 ± 0.34) (GeV/c)" 1 



(15) 



and the effective range part of the amplitude is cut off at 
TO/fTT = 1.8 GeV/c 2 . With our final fit model we have de- 
termined the preferred values of the a and r parameters 
and obtain results consistent with those given above. In- 
tegrating separately the resonant part, the effective range 
part, and the coherent sum we find that the Kq°(14:30) 
resonance accounts for 81%, the effective range term 45%, 
and destructive interference between the two terms is re- 
sponsible for the excess 26%. 

The nonresonant component of the Dalitz plot is mod- 
eled with a constant complex amplitude. We use alterna- 
tive models [27j to evaluate the model dependence of our 
results. In these studies we also use the Gounaris-Sakurai 
form [58| as an alternative model for the p°(770). 

For the angular distribution terms Tj we follow the 
Zemach tensor formalism (5tl For the decay of a 

spin B meson into a spin J resonance and a spin 
bachelor particle this gives [oil ]: 



T 



7=0 



J 



= 1 



T 



J=2 



;[3(p-q) 2 -(\p\\q\) 2 



(16) 



where p is the momentum of the bachelor particle and q 
is the momentum of the resonance daughter with charge 
opposite from that of the bachelor particle, both mea- 
sured in the rest frame of the resonance. 

The complex coefficients Cj and Cj can be parametrized 
in various ways that take the possibility of direct CP vi- 
olation into account. We have investigated the choices 
used in previous studies [13, H3] and found that they 
are susceptible to biases on the fitted parameters, par- 
ticularly when a resonant contribution is small in magni- 
tude. To ensure the good statistical behavior of our fit, 
we parametrize the complex coefficients in the following 
way: 



(xj + Axj) +i{yj + Ay 3 ) 
(xj - Axj) + i(yj — Ayj) . 



(17) 
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We have verified that these parameters have approxi- 
mately Gaussian behavior independent of their true val- 
ues. In this approach, Axj and Ayj are CP-violating 
parameters. 

To allow comparison among experiments we present 
also fit fractions (FF), defined as the integral of a single 
decay amplitude squared divided by the coherent matrix 
element squared for the complete Dalitz plot: 



FFj 



(i 



c 3 F 3 \ 



c,F,| \dm 2 K ^dm1 



\A\ +\A\ )dm 2 K7T dm: 



v 2 

7T7T 



(18) 



The sum of all the fit fractions is not necessarily unity 
due to the potential presence of net constructive or de- 
structive interference. The CP asymmetry for a given 
intermediate state is easily determined from the fitted 
parameters 



.4 



\ci\ 2 



CP, 



3 I T 1^1 

2 (x j Ax j + yjAyj) 

2 
3 



Axj+y 2 



Ay 



(19) 



The signal Dalitz-plot probability density function 
(PDF) is formed from the total amplitude as follows: 



£sig(TOLr> m L><7B) = 



(20) 



^\A\ 2 e 



^\A\ 2 e 



(\A\ 2 s+\A\ 2 e)dm 2 Kn dm^ 



where qs is the charge of the B meson candidate and 
e = e(m^,m^) and e = e(m 2 Kll ,rn 2 ^) are the signal 
reconstruction efficiencies for B + and B~ events, respec- 
tively, defined for all points in the Dalitz plot. 



III. CANDIDATE SELECTION 

The B candidates are reconstructed from events that 
have four or more charged tracks. Each track is required 
to be well measured and to originate from the beam spot. 
The B candidates are formed from three-charged-track 
combinations and particle identification criteria are ap- 
plied to reject electrons and to separate kaons and pi- 
ons. In our final state, the average selection efficiency 
for kaons that have passed the tracking requirements is 
about 80% including geometrical acceptance, while the 
average misidentification probability of pions as kaons is 
about 2%. 

Two kinematic variables are used to identify signal B 
decays. The first variable is 



AE = Eg - 



(21) 



the difference between the reconstructed CM energy of 
the B- meson candidate and y/s/2, where y/s is the total 



CM energy. The second is the energy-substituted mass 



WES = \ (s/2+p t ■ Pb) 2 /E 2 - \p B \ 



(22) 



where ps is the B momentum and (Ei,pi) is the four mo- 
mentum of the initial state, all measured in the labora- 
tory frame. The toes distribution for signal events peaks 
near the B mass with a resolution of around 2.4MeV/c 2 , 
while the AE distribution peaks near zero with a resolu- 
tion of around 19 MeV. The resolution of AE is strongly 
dependent on the position in the Dalitz plot and so in- 
stead of AE we use 



AE' 



AE 

CAB 



(23) 



where ctab is the error on AE, determined separately 
for each event. This variable exhibits no such depen- 
dence. We initially require events to lie in the re- 
gion formed by the following selection criteria: 5.200 < 
toes < 5.286 GeV/c 2 and -4.0 < AE' < 15.0. The re- 
gion of AE' below —4 is heavily contaminated by four- 
body B backgrounds and so is not useful for studying 
the continuum background. The selected region is then 
subdivided into three areas: the "left sideband" region 
(5.20 < toes < 5.26GeV/c 2 and -4.0 < AE' < 4.0) 
used to study the background AE' and Dalitz-plot dis- 
tributions; the "upper sideband" region (5.200 < Toes < 
5.286 GeV/c 2 and 7.0 < AE' < 15.0) used to study the 
background toes distributions; and the "signal region" 
(5.272 < toes < 5.286 GeV/c 2 and -4.0 < AE' < 4.0) 
where the final fit to data is performed. These three re- 
gions are illustrated in Fig. [T] Following the calculation 
of these kinematic variables, each of the B candidates is 
refitted with its mass constrained to the world average 
value of the B-meson mass Q in order to improve the 
Dalitz-plot position resolution. 



< 



15- 



10- 



5- 




5.2 5.21 5.22 5.23 5.24 5.25 5.26 5.27 5.28 5.29 

m ES (GeV/c 2 ) 

FIG. 1: Regions of the AE'-itles plane. The red/horizontal 
hatching is the "upper sideband," the green/crossed hatching 
is the "left sideband," and the blue/vertical hatching is the 
"signal region." 



The dominant source of background comes from light 
quark and charm continuum production (e + e~ — > qq, 
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where q — u,d,s,c). This background is suppressed by 
requirements on event-shape variables calculated in the 
CM frame. We compute a neural network (NN) from the 
following five variables: the ratio of the Legendre polyno- 
mial moments Lq and L2 [631 ] , the absolute value of the 
cosine of the angle between the direction of the B and the 
detector axis, the absolute value of the cosine of the an- 
gle between the B thrust axis and the detector axis, the 
output of a multivariate -B-flavor tagging algorithm [64| 
multiplied by the charge of the B candidate, and the sig- 
nificance of the measured proper time difference of the B 
decay vertices. The NN is trained using samples of off- 
resonance data and signal Monte Carlo (MC) simulated 
events. The selection requirement placed on the NN out- 
put, optimized with MC events, accepts 61% of signal 
events while rejecting 97% of background events. We see 
no significant difference between the selection efficiencies 
of positive and negative candidates. 

The reconstruction efficiency distribution over the 
Dalitz plot is modeled using two-dimensional histograms 
formed from a sample of around 24 x 10 6 i? ± — > K^n^ir^ 
MC events. All selection criteria are applied except for 
those corresponding to the invariant mass veto regions 
described below. The ratio is taken of two histograms, 
the denominator containing the true Dalitz-plot distri- 
bution of all generated MC events and the numerator 
containing the reconstructed MC events. The recon- 
structed events are weighted in order to correct for dif- 
ferences between MC and data in the particle identifi- 
cation and tracking efficiencies. In order to give better 
resolution near the edges of the Dalitz plot, where most 
reconstructed events lie, the histograms are formed in the 
"square Dalitz plot" [65[ coordinates. Linear interpola- 
tion is also applied between bins. The efficiency shows 
very little variation across the majority of the Dalitz plot 
but decreases towards the corners where one of the par- 
ticles has low momentum. The effect of experimental 
resolution on the signal model is neglected since the res- 
onances under consideration are sufficiently broad. The 
average reconstruction efficiency for events in the signal 
box for the nonresonant MC sample is 21.2%. As shown 
in Eq. (f2"0"l) , in the likelihood fit we use event-by-event ef- 
ficiencies that depend on the Dalitz-plot position and are 
different, but consistent, for B + and B~ candidates. The 
fraction of misreconstructed signal events is very small, 
~ 2%, and so such events are not treated explicitly in the 
signal model. 



IV. BACKGROUNDS 

In addition to the continuum (qq) background there 
are also backgrounds from BB events. There are four 
main sources: (i) combinatorial background from three 
unrelated tracks; (ii) three- and four-body B decays in- 
volving an intermediate D meson; (iii) charmless two- 
and four-body decays with an extra or missing parti- 
cle, and (iv) three-body decays with one or more par- 



ticles misidentified. We veto candidates from charm 
and charmonium decays with large branching fractions 
by rejecting events that have invariant masses (in units 
of GeV/c 2 ) in the ranges: 1.756 < m K + 7r - < 1.931, 
1.660 < m^+ n - < 1.800, 3.019 < m n +„- < 3.179, 
and 3.627 < m^^- < 3.747. These ranges reject de- 
cays from D° — » K + it~ (or 7r + 7r~), J/tp — > and 
ip{2S) — > ( + Cr , respectively, where £ is a lepton that has 
been misidentified. 

We study the remaining charm backgrounds that es- 
cape the vetoes and the backgrounds from charmless B 
decays with a large sample of MC-simulated BB decays 
equivalent to approximately 3 times the integrated lu- 
minosity of the data sample. Higher-statistics, exclusive 
MC samples are used to further study 68 B-meson decay 
modes and to determine the rriES) AE' , and Dalitz-plot 
distributions that are used in the likelihood fit. These 
distributions are normalized to the number of predicted 
events in the final data sample, which we estimate using 
the reconstruction efficiencies determined from the MC, 
the number of BB pairs in our data sample, and the 
branching fractions listed by the Particle Data Group [tJ 
and the Heavy Flavor Averaging Group Q. We fur- 
ther combine modes that have a similar behavior in 
each of the discriminating variables mss and AE 1 into 
a B-background category. For each category combined 
Dalitz-plot, wes, and AE' PDFs are created, and each 
is included as a separate component in the fit. The pre- 
dicted yields of BB background events in the signal re- 
gion are 619±17 (659±18) for the negatively (positively) 
charged sample. 

The background Dalitz-plot distributions are included 
in the likelihood fit through the use of two-dimensional 
histograms. For backgrounds from B decays these his- 
tograms are formed from the various MC samples. For 
the continuum background the "left sideband" data sam- 
ple is used. This data sideband also contains events 
from B decays and so MC samples are again used to 
subtract these events. Like the reconstruction efficiency 
histograms, those for the backgrounds are formed in the 
square Dalitz-plot coordinates and have linear interpo- 
lation applied between bins. Separate histograms are 
constructed for B + and B~ events. The qq and B- 
background PDFs are identical in their construction and 
the qq PDF is shown here as an example: 

Afl(m^,m^) = -(1 - q B A q q) x (24) 
/ 

, _ \ 

where A q q parametrizes possible charge asymmetry in 
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the background, and Q{m 2 Kw ,rn^„) and Q{m 2 Kv ,m^) 
are the Dalitz-plot distributions of the qq events for B + 
and B~ events, respectively. 



V. MAXIMUM LIKELIHOOD FIT 

To provide further discrimination between the signal 
and background hypotheses in the likelihood fit we in- 
clude PDFs for the kinematic variables tties and AE' , 
which multiply that of the Dalitz plot. The signal is mod- 
eled with a double Gaussian function in both cases. The 
parameters of these functions are obtained from a sample 
of K^'k^it^ MC events and are fixed in the fit to data. 
The qq tties distribution is modeled with the experimen- 
tally motivated ARGUS function 66]. The end point for 
this ARGUS function is fixed to \J~sj2 and the parameter 
describing the shape is fixed to the value determined from 
the "upper sideband." For AE' the continuum is mod- 
eled with a linear function, the slope of which is allowed 
to float in the fit to data. The BB background distri- 
butions are modeled using histograms obtained from the 
mixture of BB MC samples and are fixed in the fit. The 
yields of signal and qq events are allowed to float in the 
final fit to the data while the yield of BB background 
events is fixed. 

The complete likelihood function is given by: 



C = cxp f-J^JV*^ 



(25) 



n 



J2 NkV l k (m 2 K7T ,ml„, m ES , AE' , q B ) 



where Nk is the event yield for species k, N e is the total 
number of events in the data sample, and VI is the PDF 
for species k for event i, which consists of a product of the 
Dalitz-plot, mEs, and AE PDFs. The function — ln£ is 
used in the unbinned fit to the data. 

We determine a nominal signal Dalitz- plot model us- 
ing information from previous studies [26ll27l. l37l . l67l . [68| 
and the change in the fit likelihood value observed when 
omitting or adding resonances. In our previous study 
of B ± — * K^-k^tt^ (2(| we used a nominal model con- 
taining a phase-space nonresonant component and five 
intermediate resonant states: K*°(892)ir + , (i^7r)o°7r+, 
p°{770)K+, f (980)K+, and Xc K + . With the higher 
statistics and improved techniques of this analysis, we 
find it necessary to include additional contributions from 
^ 2 *°(1430)7r+, / 2 (1270)if+ and / x (1300)^+ in order to 
achieve a reasonable agreement of the fit with the data. 

The first of these additions improves the agreement in 
the K^ttT invariant mass projection; although some dis- 
crepancy remains, this cannot be reduced by including 
other known resonances nor by using alternative forms 
for the (Ktt)q° shape. In the tt+tt T invariant mass pro- 
jection, we find the best agreement to data is achieved by 
including both the / 2 (1270)X+ and f x (1300)K+ terms 



with X = (i.e. /x(1300) being a scalar). A reason- 
able fit can be achieved including instead a single broad 
vector resonance. However, in the case X = 1 it is nat- 
ural to identify the / x (1300) as the /9°(1450), and the 
observed large ratio of product branching fractions be- 
tweenB(B+ -> p° (U50)K+)x B(p°(U50) -> tt+tt") and 
B(B+ -» p°(770)K+) x B(p° (770) -> tt+tt-) leads us to 
conclude that this cannot be the correct physical inter- 
pretation of the data. Note that the /2(1270)if+ contri- 
bution was observed by Belle [27l ]. 

In addition, we include a small contribution 
from uj(782)K + , which is known to be present 
based on the well-measured branching fractions of 
B+ -> uj(782)K+ [H,^] and w(782) tt+tt" [zj|. Al- 
though the magnitude of this contribution is known to 
be very small, due to the narrow width of the cj(782), 
it can have a noticeable effect on the p°(770) line 
shape. Thus, our nominal signal Dalitz-plot model com- 
prises a phase-space nonresonant component and nine 
intermediate resonance states: K*°(892)ir + , (Ktt)q°tt + , 
jq°(1430)Tr+, p°(770)K+, oj(782)K+, f (980)K+, 
/ 2 (1270)JsT+, f x {1300)K+, and Y £ pif+. This model dif- 
fers from that used by Belle [27| by the inclusion of 
_fC2°(1430)7r + and the different parametrization of the 
(Kir)Q°n + and nonresonant terms. 

Our choice of parametrization of the (i^Tr)Q°Tr + and 
nonresonant terms is motivated by a number of physi- 
cal considerations. Different parametrizations result in 
changes in likelihood that are highly dependent on the 
composition of the rest of the decay model. Compared 
to our final model, including all resonant terms, the 
parametrization used by Belle gives a larger fit likelihood. 
However, this model exhibits large interference between 
the -?sTq (1430) and nonresonant terms that can be either 
constructive or destructive, leading to alternative solu- 
tions with similar likelihoods but very different values of 
the fit fractions for these terms. In the case that these 
terms interfere destructively, the sum of fit fractions far 
exceeds 100%, which we consider unlikely to be the cor- 
rect physical description. Moreover, the preferred solu- 
tion can change depending on the exact composition of 
the rest of the model - in particular, we find that this 
problem, which was previously reported by Belle (37j . is 
exacerbated by the presence of the A"2°(1430). We have 
also tried a recent proposal for the nonresonant distri- 
bution [72j that results in a worse likelihood. Therefore, 
we use in our nominal model the LASS parametrization 
for the (Ktt)q°tt + and a nonresonant term as described 
above. 

The mass and width of the /x(1300) are determined to 
be m /x = (1479±8) MeV/c 2 , T fx = (80±19) MeV, with a 
correlation of (—45 ±3)%, where the errors are statistical 
only and are determined from a fit to the 2D likelihood 
profile. These parameters are consistent with the values 
obtained by Belle m fx = (1449 ± 13) MeV/c 2 , T /x = 
(126 ± 25) MeV (statistical errors only), and with those 
listed for the / (1500) Q. 

We fit 12 753 signal candidates selected from the data 
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FIG. 2: Background-subtracted Dalitz plot of the combined 
B — > K ■k^'K data sample in the signal region. The plot 
shows bins with greater than zero entries. The area of the 
boxes is proportional to the number of entries. The depleted 
horizontal bands are the charmonium vetoes. 



using our nominal Dalitz-plot model to obtain the cen- 
tral values of the Xj, Axj, yj, Ayj parameters for each 
component. The x and y parameters of the K*°(892) are 
fixed (to one and zero, respectively) as the reference. The 
Ax and Ay parameters of the w(782) and nonresonant 
components are fixed to zero in order to improve the fit 
stability. They are allowed to vary as a cross-check and 
are found to be consistent with zero. The signal yield, qq 
background yield, and asymmetry and the slope of the 
qq AE' PDF are also floated parameters in the fit. We 
then generate a large number of MC experiments using 
the fitted values, and from the spread of results of fits 
to those experiments determine the statistical uncertain- 
ties on those parameters, as well as the central values 
and statistical uncertainties on the extracted parameters 
FFj and Acpj. This procedure takes into account cor- 
relations between the xj, Axj, yj, Ayj parameters. In 
order to make comparisons with previous measurements 
and predictions from factorization models we multiply 
each fit fraction by the total branching fraction to calcu- 
late the product branching fraction of each mode. The 
results are shown in Tables Q] and [TTJ In order to deter- 
mine the statistical significance of the direct CP viola- 
tion exhibited by a component we evaluate the difference 
A In £ between the negative log-likelihood of the nominal 
fit and that of a fit where the Ax and Ay parameters for 
the given component are fixed to zero. This is then used 
to evaluate a p- value: 



P 



2Aln£ 



f(z;n d ) dz. 



(26) 



where f(z; nj) is the x 2 PDF and rid is the number of de- 
grees of freedom, two in this case. We then determine the 
equivalent one-dimensional significance from this p- value. 
Note that this differs from the significance of Acpj ^ 0, 



since direct CP violation can be observed in a Dalitz- 
plot analysis not only through B and B amplitudes be- 
ing different in magnitude, but also by differences in their 
phases. The significance estimations were cross-checked 
using MC simulations based on the fit results. 

The if ± 7r T 7r ± signal yield is found to be 4585 ± 
90 ± 297 ± 63 and the total charge asymmetry to be 
(2.8 ± 2.0 ± 2.0 ± 1. 2)%, where the uncertainties are sta- 
tistical, systematic, and model-dependent, respectively. 
The continuum background yield and charge asymme- 
try are found to be 6830 ± 110 and (-2.8 ± 1.5)%, re- 
spectively, where the uncertainty is statistical only. The 
Dalitz plot of the data in the signal region, after subtrac- 
tion of the two background distributions (using the s Vlot 
technique [73]), can be seen in Fig. [2j Projections of the 
data, with the fit result overlaid, onto K and 7r ± 7r =F 
invariant mass distributions can be seen in Fig. [3J For 
the mx±7rT plots the requirement is made that m 1T ±^ 
is greater than 2 GeV/c 2 and vice versa in order to better 
illustrate the structures present. The agreement between 
the fit result and the data is generally very good, al- 
though the discrepancy discussed above is visible in the 
m K ± 7T T spectrum. Using the fitted signal distribution 
we calculate the average reconstruction efficiency for our 
signal sample to be 22.5%. 



VI. SYSTEMATIC UNCERTAINTIES 

The systematic uncertainties that affect the measure- 
ment of the fit fractions, phases, and event yields are as 
follows. The fixed £>i?-background yields and asymme- 
tries of the largest categories are allowed to float and the 
variation of the other fitted parameters is taken as the 
uncertainty. The effect of the statistics of the MC and 
data sideband samples used to obtain the fixed shapes 
of the efficiency, qq- and i?£?-background Dalitz-plot his- 
tograms and the SB-background histograms for toes and 
AE is accounted for by fluctuating independently the 
histogram bin contents in accordance with their errors 
and repeating the nominal fit. The uncertainties on how 
well the samples model these distributions are also taken 
into account through various cross-checks, including vari- 
ation of the charm veto range and comparison of contin- 
uum shapes between sideband and signal region in MC 
samples. 

The fixed parameters of the signal toes an d AE' 
PDFs are studied in the control sample B + — > D°ir + ; 
D° — > K + tt~. The parameters are determined from MC 
and data samples, from which biases and scale factors 
are calculated and used to adjust the parameters for the 
nominal fit. The parameters are then varied in accor- 
dance with the error on these biases and scale factors 
and the fit repeated. The uncertainties due to fixing the 
ARGUS parameter of the qq-background toes PDF are 
determined by comparing the results of the fit when the 
parameter value is obtained from off-peak data. 

To confirm the fitting procedure, 500 MC experiments 
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FIG. 3: Invariant mass projections for the data in the signal region and the fit results. The left-hand plot shows the m K ± lT 
spectrum up to 1.7GeV/c 2 . The right-hand plot shows the m 7r ± n ^ spectrum up to 1.6GeV/c 2 . The bottom plot shows the 
m 7T ± 7T ^ spectrum in the region of the XcO- The data are the black points with statistical error bars, the lower solid (red/dark) 
histogram is the qq component, the middle solid (green/light) histogram is the BB background contribution, while the upper 
blue curve shows the total fit result. For the m K ±^ plots the requirement is made that m 7r ± JT ^ is greater than 2GeV/c 2 and 
vice versa. 



were performed in which the events arc generated from 
the PDFs used in the fit to data. Small fit biases are 
observed for some of the fit parameters and are included 
in the systematic uncertainties. The contributions due 
to particle identification, tracking efficiency corrections, 
and the calculation of N bB are 4.2%, 2.4%, and 1.1%, 
respectively. The efficiency correction due to the selec- 
tion requirement on the NN has also been calculated from 
B+ —> D°tt + , D° -> K + w~ data and MC samples, and is 
found to be 0.979 ± 0.015. The error on this correction is 
incorporated into the branching fraction systematic un- 
certainties. Measured CP asymmetries could be affected 
by detector charge bias. In previous studies 0, [75[ this 
effect has been estimated to be very small compared with 
the precision of our measurements; we take it to be 0.5%. 

In addition to the above systematic uncertainties we 
also estimate effects due to model-dependence, i.e. that 
characterize the uncertainty on the results due to ele- 
ments of the signal Dalitz-plot model. The first of these 



elements consists of the parameters of the various com- 
ponents of the signal model - the masses and widths of 
all resonances, the effective range and scattering length of 
the LASS model of the (Ktt)q , the coupling constants of 
the /o(980) Flatte parametrization, and the value of the 
Blatt-Weisskopf barrier radius. The associated uncer- 
tainties are evaluated by adjusting the parameters within 
their experimental errors and refitting. The second el- 
ement is due to the different possible models both for 
the nonresonant component, which is evaluated by refit- 
ting with the parametrization used by Belle [27j | , and for 
the p°(770), which is determined by refitting with the 
Gounaris-Sakurai form. The third element is the uncer- 
tainty due to the composition of the signal model. It 
reflects observed changes in the parameters of the com- 
ponents when the data are fitted with one of the smaller 
components removed from the model and when the state 
if*°(1680) is added to the model. The uncertainties from 
each of these elements are added in quadrature to obtain 
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TABLE I: Results of fits to data, with statistical, systematic, and model-dependent uncertainties. 



Resonance 


X 


if*°(892)7r+ 


1.0 fixed 




1.718±0.084±0.064 


p°(770)K+ 


0.683±0.075±0.045 


fo(980)K+ 


-0.220±0.200±0.203 


XcoK+ 


-0.263±0.044±0.016 


Nonresonant 


-0.594±0.070±0.170 


if 2 *°(1430)7r + 


-0.301±0.060±0.030 


lo(782)K+ 


-0.058±0.067±0.018 


/a (1270) K + 


-0.193±0.043±0.022 


fx(1300)K+ 


-0.290±0.047±0.064 



Ax 



A.y 



.350 
055 

015 
073 

500 
095 

030 
014 

112 
035 

012 
134 

053 
011 

020 
033 

047 
031 



0.0 fixed 
-0.727±0.108±0.080 
-0.025±0.135±0.071 
1.203±0.085±0.052 
0.180±0.052±0.034 
0.068±0.132±0.154 
0.424±0.060±0.045 
0.100±0.051±0.010 
0.110±0.050±0.034 
-0.136±0.085±0.098 



+0.331 
-0.111 

+0.015 
-0.073 

+0.113 
-0.045 

+0.225 
-0.022 

+0.112 
-0.099 

+0.012 
-0.134 

+0.033 
-0.032 

+0.078 
-0.073 

+0.102 
-0.031 



-0.017±0.029±0.005 
-0.154±0.131±0.030 
-0.160±0.049±0.024 
-0.109±0.143±0.087 
-0.033±0.049±0.012 
0.0 fixed 

0.032±0.078±0.024 

0.0 fixed 
-0.089±0.046±0.019 

0.024±0.040±0.019 



+ 0.004 
-0.006 

+0.095 
-0.010 

+0.094 
-0.013 

+0.037 
-0.103 

+0.017 
-0.012 



+0.057 
-0.050 



+0.034 
-0.014 



+0.023 
-0.018 



-0.238±0.228±0.062 
-0.285±0.337±0.091 
0.169±0.096±0.057 
0.047±0.045±0.012 
-0.007±0.057±0.019 

0.0 fixed 
0.007±0.086±0.017 

0.0 fixed 
0.125±0.058±0.021 
0.056±0.087±0.044 



+ 0.144 
-0.018 

+0.221 
-0.019 

+0.133 
-0.027 

+0.046 
-0.018 

+0.006 
-0.111 



+0.025 
-0.034 



+0.034 
-0.025 

+0.010 
-0.036 



the final model-dependence. 

VII. SUMMARY AND DISCUSSION 

Our results are shown in Tables U and [TTJ The to- 
tal branching fraction B(B + — > K + tt~tt + ) = (54.4 ± 
1.1 ± 4.5 ± 0.7) x 10"~ 6 is compatible with Belle's mea- 
surement of (48.8 ± 1.1 ± 3.6) x 10" 6 [13. This result 
was cross-checked by using the same procedure to mea- 
sure the B + —>■ D°ir + ; D° — > K + ir~ branching fraction, 
which was found to be consistent with the PDG value Q • 
The total charge asymmetry for B + — ► K + tt^-k + has 
been measured to be consistent with zero to a higher 
degree of precision than previous measurements. 

We see evidence of large direct CP violation in 
B + — > p°(770)K + , consistent with the findings of our 
previous analysis [26| and of Belle [27| . The statistical 
significance of the direct CP violation effect is found to 
be 3.7<7 from the change in likelihood when the Ax and 
Ay terms associated with p°(770)K + are fixed to zero. 
We have verified this estimate of the significance using 
MC simulations. As experimental systematic uncertain- 
ties are much smaller than the statistical errors, they do 
not affect this conclusion. We have cross-checked the ef- 
fect of the choice of the Dalitz model on the significance. 
We find that the significance remains above 3<r with al- 
ternative models, including that used by Belle [27| . 

Plots of the 7r + 7r~ invariant mass projections in the re- 
gion of the p°(770) and / (980) are shown separately for 
B + and B~ candidates in Fig. 0J In an attempt to high- 
light the direct CP violation effect, we also show plots 
with the data further subdivided into positive and nega- 
tive values of cosOh — P ■ q / (\p\W\), using the notation 
of Eq. (fTC)) . The asymmetry in the excess of events above 
background in the p°(770) region is particularly apparent 
in the distributions with the requirement cos9h > 0. 

The statistical significance of direct CP violation in 
B + — > /2(1270)AT + is also above 3c, but this result suf- 
fers from large model uncertainties. The AT*°(892)7r + , 



(Ktt)q°tt + , and AT2°(1430)7r + charge asymmetries are all 
consistent with zero, as expected. 




0.6 0.7 0.8 0.9 1 1.1 0.7 0.8 0.9 1 1.1 1.2 



(GeV/c 2 ) m m (GeV/c 2 ) 

FIG. 4: Projection plots of the 7r + 7r _ invariant mass in the 
region of the p°(770) and /o (980) resonances. The left (right) 
plots are for B~ {B + ) candidates. The top row shows all 
candidates, the middle row shows those where cos 9h > 0, 
and the bottom row shows those where cos 9h < 0. The data 
are the black points with statistical error bars, the lower solid 
(red/dark) histogram is the qq component, the middle solid 
(green/light) histogram is the BB background contribution, 
while the blue open histogram shows the total fit result. 



After correcting for the secondary branch- 
ing fraction B(K*° (892) K+n-) = § we find 

B{B+ -> K*°(892)tt+) to be (10.8±0.6± 1.1 tot) x 10 ~ 6 - 
Similarly we find B(B+ -> / 2 (1270)A'+) to be (0.88 ± 
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TABLE II: Summary of measurements of branching fractions (averaged over charge conjugate states) and CP asymmetries. 
Note that these results are not corrected for secondary branching fractions. The first uncertainty is statistical, the second 
is systematic, and the third represents the model dependence. The final column is the statistical significance of direct CP 
violation determined as described in the text. 



Mode 




Fit fraction (%) 


B(B+ -> Mode)(10" 6 ) 


Acp (%) 


DCPV sig. 


K + n~n + total 








54.4 ± 1.1 ±4.5 ±0.7 


2.8 ± 2.0 ± 2.0 ± 1.2 




A"* (892)7r+; A"*°(892) -> 




13.3 ±0.7 ±0.7: 


h0.4 
-0 9 


7.2 ±0.4 ±0.7^;| 


±3.2 ± 5.2 ± 1.1 to% 


0.9cr 


(/f7r)3°7r + ; (#7r)8° -> K+ 




45.0 ± 1.4 ± 1.2 1 


-12.9 
-0 2 


24.5 ± 0.9 ± 2.1 


±3.2 ± 3.5 ± 2.0 


1.2(7 


p°(770)A+; p°(770) -> rr 4 




6.54 ±0.81 ±0.58 


+0.69 
q 26 


3.56 ±0.45 ±0.43^;^ 


+44 ± 10 ±4 If 3 


3.7cr 


/o(980)A+; / (980) rr 4 




18.9 ± 0.9 ± 1.71 


h2.8 
-0.6 


10.3 ± 0.5 ± 1.3 to.l 


-10.6 ± 5.0 ± lTtto 


1.8(7 


XcoK + ; Xco -> Tr+TT" 




1.29 ±0.19 ±0.15 


+0.12 
-0.03 


0.70 ± 0.10 ± 0.10 ±g;g| 


-U±15±3±l 


0.5cr 


X + 7r _ 7r + nonresonant 




4.5 ± 0.9 ± 2.4 1 


0.6 
1.5 


2.4 ± 0.5 ± 1.3 ±g;| 






ifl°(1430)7r + ; A- 2 *°(1430) 




3.40 ± 0.75 ± 0.42 


+0.99 
-0.13 


1.85 ±0.41 ±0.28^;^ 


+5 ±23 ±4 


0.2cr 


u)(7$2)K + ; o;(782) -> Tr+Tr" 


0.17 ±0.24 ±0.03 


+0.05 
-0.08 


0.09 ±0.13 ±0.02 t l 






/ 2 (1270)^+; / 2 (1270) -> 


7T + 7T 


0.91 ±0.27 ±0.11 


+0.24 
-0.17 


0.50 ±0.15 ±0.07^;^ 


-85 ±22 ±13 If 


3.5cr 


/ X (1300)JC+; /x(1300) -> 


7T + 7T~ 


1.33 ±0.38 ±0.86 


+0.04 
-0.14 


0.73 ± 0.21 ± 0.47 1 ° 2 8 


+28±26± 13 ±1 


0.6cr 



0.26±0.13 to°ol) x 10- 6 and B(B+ -» if 2 *°(1430)7r+) 
to be (5.6 ± 1.2 ± 0.8 1\% ± 0.1) x 10" 6 , where the fourth 
errors are due to the uncertainties on the secondary 
branching fractions. These and the other branching 
fraction measurements are, in general, consistent with 
previous measurements. The product branching fraction 
of B+ -> ui(782)K+; w(782) -> tt+tt" agrees with the 
expectation from previous measurements, albeit with 
large uncertainties. The only disparities between our 
results and those reported by Belle [27| arise from 
the different treatments of the K + ir~ S-wave. The 
forward-backward asymmetry apparent in both the 
K*°(892) and / (980) bands in Fig. lis well reproduced 
by the fit and is due to S-P-wave interference in the 
Dalitz plot. 

The (Ktt)q component is modeled in our analysis by 
the LASS parametrization [56| . which consists of a non- 
resonant effective range term plus a relativistic Breit- 
Wigner term for the Kq°(1430) resonance itself. This 
parametrization makes use of the available experimental 
information, however the size of the phase space here 
means that the range of K + n~ invariant masses un- 
der consideration is much greater than in previous stud- 
ies. The agreement between the model and the data in 
the region of the A"q°(1430) is not as good as in the 
remainder of the Dalitz plot but all alternative mod- 
els tried yielded poorer results. If we assume that the 
model used is correct then we can calculate the branch- 
ing fraction for B + -> K$°(U30)tt + and find it to be: 

(32.0 ± 1.2 ± 2.7 ti l ± 5 - 2 ) x icr6 > where the fourth 
error is due to the uncertainty on the branching frac- 
tion of Kq°(1A3Q) — > Ktt combined with the uncertainty 
on the proportion of the (Ktt)q° component due to the 
Kq°(1430) resonance. In addition we can calculate the 



total nonresonant contribution by combining coherently 
the nonresonant part of the LASS parametrization and 
the phase-space nonresonant. We find this total nonreso- 
nant branching fraction to be: (9.3 ± 1.0± 1.2 ±££ ± 1.2) x 
10~ 6 , where the fourth error is due to the uncertainty on 
the proportion of the (i^7r)o° component that is nonres- 
onant. The Belle collaboration finds somewhat larger 
i^Q°(1430) and nonresonant branching fractions though 
they treat the A~o (1430) as a Breit-Wigner component, 
separate from the rest of the S-wave modeled as a non- 
resonant amplitude that has variation in mag nitude but 
no variation in phase over the Dalitz plot [27] . 

In conclusion, we have performed a Dalitz-plot anal- 
ysis of B ± — > K ± tt t tt ± decays based on a 347.5 fb _1 
data sample containing (383.2 ± 4.2) x 10 6 BB pairs col- 
lected by the BABAR detector. To obtain a good fit to 
the data we find that contributions from /2(1270)A: + 
and /x(1300)A' + are necessary, with /x(1300) being 
a scalar with parameters m/ x = (1479 ± 8)MeV/c 2 , 
r /x = (80 ± 19) MeV, with a correlation of (-45 ± 3)%, 
where the errors are statistical only. These are consis- 
tent with the mass and width of the /o(1500). In the 
A' ± 7r =F invariant mass projection some discrepancy re- 
mains in the 1200-1400 MeV/c 2 range - our model in- 
cludes (A"7r)g°7r + (using the LASS parametrization) and 
A:2°(1430)7T + terms in this region, although some alter- 
native models also give a good description of the data. 
We measure CP-averaged branching fractions and direct 
CP asymmetries for intermediate resonant and nonreso- 
nant contributions. Our results are consistent with the 
standard model, and can be used together with results 
from other B — > A:7T7t decays to obtain constraints on 
the CKM phase 7. We find evidence for direct CP viola- 
tion in the decay B^ — > p°(770)A: ± , with a CP- violation 
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parameter A C p = (+44 ± 10 ± 4lf 3 )% . These results 
supersede those in our previous publication [26| . 
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